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^ . Abstract 

<: 

We use high-spatial resolution (~ 100 pc) rotation curves of 83 spiral galaxies 
to investigate the mass distribution of their innermost kpc. We show that, in this 
region, the luminous matter completely accounts for the gravitational potential and 
no dark component is required. The derived I-band disk mass-to-light ratios yj agree 
well with those obtained from population synthesis models and correlate with color 
00 . in a similar way. We find strict upper limits of ~ 1O 7 M0 for the masses of compact 

bodies at the center of spirals, ruling out that these systems host the remnants of 
the quasar activity. 

^: 

P ■ 1 Introduction 

£3 ■ 

Ctf ■ 

There is increasing direct evidence that, at the centers of bulge-dominated 
^ . galaxies there reside Massive Dark Objects (MDO), which are probably the 

remnants of the engines that once powered the QSO phenomenon (Ho, 1998; 
Kormendy & Richstone, 1995). In fact, virtually every hot galaxy hosts a 
MDO/BH with a mass ranging from ~ 1O 8 M to 2 x 1O 1O M , similar to those 
related to the QSO phenomenon. For disk galaxies, the situation is different 
and much more uncertain. A direct determination of the central mass has been 
obtained only in very few cases which include our own Galaxy where a black 
hole with a mass 2 x 10 6 M Q has been discovered (Ghez et al.1999, Genzel 1998; 
see also Salucci et al 1999 for other few cases). Remarkably, these masses do 
not exceed 1O 7 M ; however, the lack of detections of very massive objects 
(M MDO > 1O 8 M ) cannot be ascribed to observational biases. In fact, the 
rotation curves (RC), in great number available down to r in ~ 100 pc, could 
easily expose central bodies with masses of the order of ~ 1O 8 M , given that 
usually V(ri n ) ~ 10 km/ s, or equivalently, the stellar mass inside lOOpc barely 
reaches 1O 7 M . On the other hand, the same RC analysis determining the 
MDO mass, obtain also the mass distribution of the innermost regions. This 
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is particularly important in the case of low luminosity DM-dominated spirals 
(see Persic and Salucci, 1990) for which it is generally difficult to disentangle 
the disk component from the whole bulk of gravitating matter. In any cases, let 
us stress that to proper model the region where the luminous matter dominates 
is indispensable to infer the structural properties of dark matter. 

The aim of this work is to study the innermost kinematics of spirals in order to 
a) investigate the distribution of matter in relation to that of light b) set strin- 
gent upper limits on the MDO masses. The plan of the paper is the following: 
in section 2 we analyse a sample drawn from the 967 rotation curves (RC's) 
of Persic and Salucci (1995), in section 3 we derive and discuss the disk prop- 
erties and in sections 4 and 5 we derive upper limits on the MDO/BH masses 
and discuss the results. In this paper we assume H Q = 75 km/s Mpc^ 1 , and 
take, as the reference magnitude in the B-band, M* = —20.5, which translates 
to M* — 21.9 in the I-band (e.g. Rhee, 1997). All luminosities correspond to 
the I-band unless otherwise specified. 



2 Inner Rotation Curves of Late Types Spirals and mass modeling 

Recently, about a thousand rotation curves of spiral galaxies (PS95), tracing 
the kinematics inside the central kpc, have been available. Among these, about 
one hundred have at least one measurement < 200 pc and several < 500 — 
1000 pc. For these objects the detection of a MDO of ~ 1O 8 M or more, if 
present, is guaranteed and so is the determination of the stellar disk mass, 
provided that in this region the stellar disk is the major mass component. 

Therefore, we select from PS95 the largest sample of rotation curves equally 
distributed in magnitude interval and of sufficient high-quality for our pur- 
poses, by setting the following criteria: i) each RC has at least 4 measurements 
inside a radius of 340(1 — ( M/ + 15 ) ) p C ^ (350 — 500) pc, ii) for each RC, the 
innermost data is situated at a radius r in < 100(1 — ( Ml + 15 ^ ) ~ (100 — 150) pc. 
Let us notice that to restrict the criteria reduces the number of RC's with- 
out improving their already high quality that however rapidly decreases as 
they are relaxed. The rotation curves resulting from the selection are gen- 
erally smooth, axisymmetric and with negligible non circular perturbations. 
It is evident that there is a very good agreement between the rotation fields 
on the receding and the approaching arm (see PS95). The sample (Sample B 
in Salucci et al., 2000) has 83 objects, well mixed in luminosity and almost 
equally distributed between Sb-Sc and Sd-Im Hubble types. 

We aim to reproduce the RCs by a mass model featuring: (i) a Freeman disk 
of length-scale Rd, (derived in PS95) and mass M D) which contributes to the 
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Fig. 1. Mass models of the 83 objects oj^ sample B. The dashed lines indicate the 
0D model (in a few bulge has been included) while the solid line indicates 

the disk+BH model. 



circular velocity as: 

VKR) = \GM D /R D x 2 (I K - (1) 

where I n ,K n are modified Bessel functions and x = R/Rd, and (ii) a MDO 
which contributes to V(R) as: 

Vl 1DO {R) = GM MDO /R (2a) 

where the MDO mass is: 

M MDO = fG- l V\r m )r m (26) 

with < / < 1. 

In Fig.(l) we compare the velocity data with the best-fitting mass models in 
the cases of 1) / = (only disk, OD) mass model and of 2) / > (disk+ 
black hole) mass model (in practice: / ~ 1). As result, we find no central- 
body-dominated (CBD) rotation curve, i.e. no RC shows the Keplerian fall-off 
expected for a CBD RC. On the contrary, rotation curves are strikingly close 
to those predicted by the OD model, i.e. by a self-gravitating exponential thin 
disk with radially constant mass-to-light ratio. The results of the mass mod- 
eling are: in 57 objects the OD velocity curve accounts for the rotation data 
in a excellent way and it lies within the data error bars (for R < Ribd see 
below). In 14 cases(3]the luminous matter accounts for the rotation curve in a 
satisfactory way: data and model differ utmost by 2a and no evidence for an 
additional component emerges. In 9 cases' the OD fits are reasonable expe- 
cially if we consider that the corresponding RC's have non-negligible internal 
dispersion and/or some asymmetry. The low-luminosity galaxy 545- G3 is dom- 
inated by dark matter at any radius. Finally, in 2 cases; 346-G26, 410-G19, 
the OD model reproduces the data with some difficulty 

A simple inspection of the results shows that in every galaxy only for R > 
-RfBnQthe dark component begins to significantly contribute to the total grav- 
itational potential. We then confirm the picture of Salucci & Persic (1999) 
and Salucci etal, (2000) according to which, the luminous matter dominates 
an innermost region of spirals of size Ribd- This truly baryonic scale turns 
out to be a function of galaxy luminosity by ranging from 0.5Rd ~ 0.5kpc 
to 2Rd ~ 30kpc along the luminosity sequence Such a transition can be 
seen, e.g. in the RC of 55-G4: OD curve reproduces the data only out to 
R = Ribd ~ 0.5kpc. More in general, the best-fit models (see Fig 1) clearly 

1 349-G6, 38-G23, 35-G18, 157-G20, 436-G3, 468-G23, 547-G24, 550-G7, N4348, 
358G-63, 582-G12, 547-G31, 196099, N4705 

2 347-G84, 328-G43, 249-G16, 84-G10, N755, 79-G14, M222025, N7218, 249-G16 

3 Inner Baryon Dominance 
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Fig. 2. The 3D space of galaxy properties. Md — {B — I) — Lj. 

show that, with the exception of the least luminous galaxy, at R ~ 0.5 kpc, the 
DM contribution to the RC lies below the detection threshold. On the other 
hand, in order to become the major component at a 2 — 3 disk scale-lenghts, 
the dark mass, negligible at 0.5kpc must strongly increase with radius. 

We conclude by stating that the OD mass models reproduce the innermost 
kpc of the rotation curves extremely well: on this scale, there is it no hint of 
a dark component. 



3 Disk properties 

Disk masses are derived from eq (1); Although the uncertainties on the disk 
mass-to-light ratio values must include model-fit uncertainties, photometry 
errors, and uncertainties on the assumed inclination and distance (the lat- 
ter present also for Tully-Fisher distances). The uncertainty budget can be 
roughly estimated as: 5M D /Lj ~ 0AM D /Lj. 

The present sample of DM-free mass models is a most suitable one for investi- 
gating the disk masses in relation with the photometric properties of galaxies. 
In detail, this is done with a sub-sample of Sample B comprising 28 objects 
for which /-band luminosities and B — I colors are available. In Table 1 we 
report best-fit disk masses, color, mass-to-light ratios M^/Li and luminosity 
for these objects. B-I color are good indicators of stellar mass-to-light ratio 
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Fig. 3. Logarithm of the mass-to- light ratio in the I band as a function of the B-I 
color. The solid line is from Vazdekis etal. 1996. The uncertainty in B-I is 0.1 mag 

(Vazdekis et al 1996). In Fig. 2 we the show yj vs. B — I colors alongside with 
with the same relation predicted by population synthesis models. We find 

3^0.63x(H^) 2 (3a) 

in good agreement with Vadzekis et al Finally, we compare the kinematical disk 
masses M D with those derived from the galaxy specto-photometry: M phot = 
yi {B — I) Lj (let us remind that yi is the average stellar mass-to-light ratio 
for a stellar population of color B — I. As a quite general outcome of the past 
history of galaxies, Vazdekis et al. 1996 found 

y T = b (B - I) + c (36) 

where b ~ 0.8 and the value of the constant c is irrelevant. We find: 

M D = (0.96 ± 0.1)M phot (4) 

(28 d.o.f., see Fig 4), i.e. dynamical and photometric mass estimates statisti- 
cally coincide. Notice that the disk mass and I-band luminosity are not directly 
proportional: 

log M D = (0.79 ± 0.1) logLi + const (5) 

with the log slope being significantly different from unity. Then, the existence 
of a relationship between the maximum rotation velocity and the galaxy lu- 
minosity is far from being a trivial one, in that the J-band luminosity is not a 
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Fig. 4. Kinematical mass Md vs photometric disk masses. Units in logM@ 

straightforward measure of the stellar mass. It is worth to note that in the 3-D 
space defined by the (logLj, log(M D /Lj), B — I) coordinate vector, spirals 
are not randomly distributed, but occupy a very thin plane (see Fig. 2). 



4 MDO/BH masses. Upper limits 

The very good spatial resolution of the RC's of our sample plays a crucial 
role in the determination of strict and reliable MDO upper limits. In fact, 
since at small radii V oc R and the mass inside r in scales as rf n , a spatial 
resolution of Ti n ~ lOOpc yields to a mass resolution of ~ lO^AfQ, given that 
V{r in ) ~ 10km/ s. In detail, the ratio of the rotational velocities (VMDo/Vd) 2 oc 
(M MDO /r in ) 3 , implies that a higher spatial resolutions lead to stronger limits 
on Mmdo ( or easier detection). 

We now force the presence of a central MDO/BH by adopting the maximum 
possible MDO mass compatible with the rotation curve. We determine the 
upper mass limits by assuming / = 1 or / = f max with / = f max the value that 
brings the model RC la higher than the innermost data. Note that for / = 1 
all the mass inside r in is in the black hole. We typically find: f max = 0.7 — 0.8 
and that, in spite of having an additional parameter, the disk+black hole 
model performs significantly worse than the OD model. 
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Notice that only in a few cases f max ~ 0.2, furthermore, no result changes if 
we assume f max = 1 for all the objects. M BH is then obtained by substituting 
in equation (2b) the corresponding values of /. These are shown in Table 1: it 
is evident that typically M MDO « 1O 8 M especially at low luminosity. It is 
clear that these limits are short of the mass residing in a accreting BH which 
is required to power high redshifts quasars, thus, spirals do not host QSO 
remnants (see also Salucci etal 2000). More specifically, if central MDO's of 
masses ~ 2 M MDO (and still < 1O 8 M ) were actually present at the center 
of the late type spirals of our sample, they would have affected the inner kpc 
of the available RC's more strongly than MDO's of ~ 1O 9 M affect the inner 
kinematics of ellipticals (Magorrian et al., 1998). 

In Fig. 5 we show the relationship M BH vs the stellar disk mass M D , and, as 
a reference, we plot the extrapolation to spirals of the elliptical's stellar mass- 
MDO/BH mass relationship M MDO = 0.005 x M stars (Magorrian et al.,1998, 
Kormendy & Richstone 1995 ). We immediately realize, that in spirals, the 
MDO masses, if different from zero, are however much smaller than those de- 
tected at the centers of ellipticals of same stellar luminosity; whether this is 
related to their less prominent bulge content or does reflect also a morpholog- 
ical mass segregation, will be considered elsewhere (Salucci et al 2000). Notice 
that in the same paper we perform a detailed investigation on possible biases 
occurring in the present estimate of the MDO upper limits. 
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5 DISCUSSION 



Accurate mass models of the innermost kpc of 83 spirals derived from high 
spatial resolution RC's, reveal that, inside this region, the luminous matter 
fully accounts for the observed kinematics: a stellar disk of constant mass-to- 
light ratio (in some cases in conjunction with a spheroid) is virtually the only 
mass component. In this region, with the exception of very high luminosity 
galaxy with complex dynamics, we find 

V 2 (R) oc light(< R)/R (6) 

which is indicative of the absence of a sizable dark component. Therefore, the 
only possibility for the presence of a dark halo in the central regions of spirals 
is that it conceals itself below the detection threshold. 

More in detail, inside the innermost kpc or so, the contribution to the grav- 
itational potential from a dark component cannot exceed 10%; this is in dis- 
agreement with claims and scenarios in which a dark halo has a major role at 
any radii. Moreover, in the same objects studied here, the DM, undetected at 
R < 1/3 Rd, begins to dominate the mass distribution at R ~ SRc (PSS). 
Then, the dark mass inside a giv en radius R must increase with R very steeply, 
e.g. as in the case of a constant density distribution. 

In a DM-free environment, we have been able to estimate the disk masses also 
for objects in which the global influence of the dark matter is relevant and to 
relate disk colors and masses across a large range of galaxy luminosities. The 
next step will be to pursue a coordinated kinematical and photometric study 
to tackle open cosmological issues. 



6 Appendix 

In the appendix we produce the two Tables of galaxy and MDO/BH properties. 
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Table 1 

(1) Galaxy name (2) Logarithm of the disk mass in solar units (3) B-I color 
Logarithm of the disk mass-to- light ratio (I-band). (5) I-Band luminosity. 
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556-g23 


7.3 


582-gl2 


7.7 


362-gll 


6.2 


563-gl4 


7.3 


79-gl4 


7.6 


374-g8 


5.2 


566-g22 


6.5 


79-g3 


6.2 


380-g23 


6.8 


575-g53 


6.0 


196099 


6.9 


406-g26 


5.2 


58-g30 


6.6 


m-222025 


8.6 


410-gl9 


7.0 


60-g25 


6.4 


n4705 


7.2 


418-gl 


5.9 


84-glO 


5.5 


n697 


8.4 


1407 


7. 


196099 


6.9 







Table 2 

BH/MDO upper limits for the galaxies of the sample, in units of logM Q 
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